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We use Monte Carlo and molecular dynamics simulation to study a magnetic tip-sample interac- 
tion. Our interest is to understand the mechanism of heat dissipation when the forces involved in the 
system are magnetic in essence. We consider a magnetic crystalline substrate composed of several 
layers interacting magnetically with a tip. The set is put thermally in equilibrium at temperature 
T by using a numerical Monte Carlo technique. By using that configuration we study its dynamical 
evolution by integrating numerically the equations of motion. Our results suggests that the heat 
dissipation in this system is closed related to the appearing of vortices in the sample. 

PACS numbers: 75.10.Hk, 75.30.Kz, 75.40.Mg 



I. INTRODUCTION 

The friction between two sliding surfaces is one of the 
oldest phenomena studied in natural sciences 1]. In a 
macroscopic scale it is known that the force of friction 
between surfaces satisfies some rules: 1 - The friction is 
independent of the contacting area between surfaces. 2 - 
It is proportional to the normal force applied and 3 - The 
force of kinetic friction is independent of relative speed 
between surfaces. That behavior being the result of many 
microscopic interactions between atoms and molecules of 
both surfaces, it is also dependent on the roughness, tem- 
perature and energy dissipation mechanisms. Therefore, 
to understand friction it is necessary to understand its 
microscopic mechanisms iJ]. 

For several applications of now a days technology the 
understanding of how heat is dissipated when mobile 
parts are involved, plays an important role. The avail- 
ability of refined experimental techniques makes it now 
possible to investigate the fundamental processes that 
contribute to the sliding friction on an atomic scale. Is- 
sues like how energy dissipates on the substrate, which 
is the main dissipation channel (electronic or phononic) 
and how the phononic sliding friction coefficient depends 
on the corrugation amplitude were addressed , and par- 
tially solved, by some groups. |3jI31- Less known is the 
effect on friction of a magnetic tip moving relative to a 
magnetic surface. Applications of sub-micron magnets 
in spintronics, quantum computing, and data storage de- 
mand a huge understanding of the sub-micron behavior 
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of magnetic materials. The construction of magnetic de- 
vices has to deal with distances of nanometers between 
the reading head and the storage device. That makes 
the study of tribological phenomena crucial to under- 
stand and produce technologically competitive devices 
Hi M, H^ • In particular the dissipation of heat in mag- 
netic dispositives is a very serious problem. For example, 
in a magnetic hard disk for data storage the reading head 
passing close to the surface of the disk transfers momen- 
tum to it. That momentum transference rises locally 
the temperature. Depending on the rate transfer and 
the capability of the disk to transfer that energy to the 
neighborhood some, or all, the information stored in the 
disk can be lost. In the last decade, the progress in the 
magnetic recording media and the reading head technol- 
ogy has made the recording density doubled almost every 
two years. The magnetic bit size used in the most ad- 
vanced hard-disk-drives is as small as 0.5 x O.Obm^, while 
a giant magneto-resistance head is used to read the bit. 
This magnetic bit size can still be diminished by using 
materials with high magnetic anisotropy. However, the 
paramagnetic limit, when the anisotropy energy becomes 
comparable to the thermal fluctuations, is not far to be 
attained. That is a physical limitation to the technology 
for the production of magnetic recording media. In that 
case, the comprehension of the microscopic mechanism 
of heat dissipation is crucial. Since friction is an out- 
of-equilibrium phenomenon its study presents numerous 
theoretical and experimental challenges. 

In magnetic films for which the exchange interaction 
is less than the separation between layers form quasi 2 
dimensional (2d) magnetic planar structures. In general 
the magnetization of such films is confined to the plane 
due to shape anisotropy. An exception to that is the 
appearing of vortices in the system. A vortex being a 
topological excitation in which the integral of line of the 
field on a closed path around the excitation core precess 



by 27r (— 27r) (See figure n). To the purpose of avoid- 
ing the high energetic cost of non-aUgned moments, the 
vortices develop a three dimensional structure by turn- 
ing out of the plane the magnetic moment components in 
the vortex core, llj For data storage purposes, magnetic 
vortices are of high interest since its study provides fun- 
damental insight in the mesoscopic magnetic structures 
of the system 12J. 

In this paper we use a combined Monte Carlo- 
Molecular Dynamics (MC-MD) simulation to study the 
energy dissipation mechanism in a prototype model con- 
sisting of a reading head moving close to a magnetic 
disk surface. A schematic view is shown in figure |21 
Our model consists of a magnetic tip (The reading head) 
which moves close to a magnetic surface (The disk sur- 
face.). The tip is simulated as a cubic arrangement 
of magnetic dipoles and the surface is represented as a 
monolayer of magnetic dipoles distributed in a square 
lattice. We suppose that the dipolc interactions arc 
shielded, so that, we do not have to consider them as 
long range interactions. This trick simplifies enormously 
the calculations putting the cpu time inside reasonable 
borders. The dipole can be represented by classical spin 
like variables S = S^x + Syy + SzZ. The total energy of 
this arrangement is a sum of exchange energy, anisotropy 
energy and the kinetic energy due to the relative move- 
ment between the tip and the surface as follows. 
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FIG. 1: Schematic view of a vortex (left) and anti- vortex 
(right) for spins of the same length. 



between the tip and the surface. The symbol < ij > 
means that the sums are to be performed over the first 
neighbors. For the tip-surface interaction, we suppose 
that the coupling, Jh-s, is ferromagnetic. By consider- 
ing that Jh-s is a function of distance, will allow us to 
study the effects of the relative tip-surface movement. 
The exchange anisotropic term A, controls the kind of 
vortex which is more stable in the lattice. There is a 
critical value of the anisotropy, Ac ~ 0.7 J |T2|, such that 
for A < Ac the spins inside the vortex core minimizes 
the vortex energy by laying in an in-plane configuration. 
For A > Ac the configuration that minimizes the vortex 
energy is for the spins close to the center of the vor- 
tex to develop a large out-of-plane component. The site 
anisotropy, Z?, controls the out-of-plane magnetization 
of the model. It is well known that a quasi-two dimen- 
sional system with interaction as in equation|21undergoes 
a BKT phase transition for Dh sufficiently small (In gen- 
eral D/ J « A.), at some critical temperature Tbkt- If, 
D/J = A = 0, Tbkt « 0.7. For D/J » A the system 
has a second order phase transition at Tc which depends 
on D. 

In the sections below we will discuss the importance of 
vortex formation for the energy dissipation of a magnetic 
material. In the section|nlwe introduce general aspects of 
the numerical method used, in the section ITITI we discuss 
our results and in section Hvl we present our conclusions. 
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Uh-s = - ^ Jh-s {\rhi - rsj\) {Sm ■ Ssj) (4) 



Jh-s {\rht - rsj\) = Jo exp{-Q: (rhi-sj - tq) } (5) 

In equation n the first term, 'pff^_s)i/2'm(h)i , stands for 
the relative kinetic energy : surface-reading head (s- 
h). The second term, Uspin , accounts for the magnetic 
dipoles interactions: in the tip {Uh) and in the surface 
(Us). The last term, Uh-s, is the interaction energy 



II. SIMULATION BACKGROUND 

In this section we describe the numerical approach we 
have used to simulate the model. The simulation is done 
by using a combined Monte Carlo-Molecular Dynamics 
(MC-MD) procedure. The particles in our simulation 
move according Newton's law of motion which can be 
obtained by using the hamiltonian ^ The spins evolve 
according to the equation of motion 13] 
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spectively. In the second set, we put Ds = Dh = and 
As = A/i = 0.6 J, 0.9 J to get vortices with in-plane and 
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FIG. 2: Schematic view of our simulation arrangement. The 
tip is simulated as a 5 x 5 rigid lattice of particles each carrying 
a spin S. The surface is simulated as an arrangement of size 
20 X 20. The tip is allowed to slid over the surface with initial 
velocity v = vox. Periodic boundary conditions are used for 
the in plane directions of the surface. The arrows represent 
the spins directions. 



As we are mainly interested in the magnetic effects we 
consider the particles as fixed letting the tip slid over the 
surface. This generates a set of 3A^-|-3 coupled equations 
of motion which are solved by increasing forward in time 
the physical state of the system in small time steps of size 
St = 10~^J~^. The resulting equations are solved by us- 
ing Runge-Kutta's method of integration jlj,[i^ [13 ■ The 
surface is arranged as a rigid 20a x 20a square lattice with 
periodic boundary conditions in the xy direction, where 
a is the lattice spacing. The head is simulated as a rigid 
5a X 5a square lattice. With no loss of generality the 
lattice spacing will be taken as a = 1 from now on. In 
figure 121 we show a schematic view of the arrangement 
used in our simulation. Initially we put the tip at a large 
distance from the surface, so that, Jh-s = 0. By using 
the MC approach we equilibrate the system at a given 
temperature, T. By controlling the energy of the system 
we consider that the system is in equilibrium after 10"" 
MC steps. Once the thermal equilibrium is reached we 
observe the evolution of the separated systems for a small 
time interval, r, for posterior comparison. After that an 
initial velocity, v = vqx is given to the tip. We follow the 
system's evolution by storing all positions, velocities and 
spin components at each time step for posterior analy- 
sis. A quantity of paramount importance for us is the 
vortex density, ipv{t), at the surface layer. We calculate 
the vortex density as a function of temperature and time 
for several values of the parameters A^, Xh and Ds, D^- 
As discussed in the following the vortex density will be 
related to energy dissipation. The energy is measured in 
units of Jh ~ Js — J 1 temperature in units of J /kB, time 
in units of J^^ and velocity in units of a/ J~^ where fcs 
is the Boltzmann constant. 



III. RESULTS 

The system is simulated for several temperatures and 
anisotropy parameters. In all of them we have fixed 
Jh = Js = J, Jo = 2 J and a = 1. In the first set of 
simulations we use A^ = A/j = IJ and the site anisotropy 
Ds = Dh = O.IJ, — O.IJ to obtain a system with out- 
of-plane (Ising-like) and in-plane (XY) symmetries re- 
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FIG. 3: Vortex Density as a function of the temperature. 
The symbols are for Xa = Xh = 1 with Dh — Ds — O.IJ (o), 
Dh = Ds = -O.IJ (D) and Ds =Dh=0 with A^ = A^ = 0.6 J 
(0) and Ah = As = 0.9J (a). The lines are only guides to the 
eyes. 

evolution, we have to know how the vortex density, ify, 
depends on the temperature. In figure 13 we show a plot 
of ifv as a function of the temperature for the simulated 
models. The vortex density increases monotonically with 
T. If there is any relation between vortices and energy 
dissipation, it is natural to think that an increase in the 
vortex density is related to an increase in temperature. 

With the head far from the surface, we start the time 
evolution of the system at i = with f = 0. This part of 
the simulation serves as a guide to the rest of the simula- 
tion. Only thermal fluctuations of the vortex density can 
be seen. At t = 200 the tip is released with initial veloc- 
ity vq. For t > 200, the reading head starts to interact 
with the surface. Some kinetic energy is transferred from 
the head to the surface and we expect the vortex density 
to increase. We will see in the following that depending 
on the initial conditions and the symmetry of the sys- 
tem (See equation^ several things can happen. In 
the figures 01 to [7| we show plots of the averaged vortex 
density as a function of the time for the set of the simu- 
lated parameters as discussed above. In each figure, the 
graphics from left to right and top to bottom correspond 
to r = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7 and 0.8. 

For the system with out-of-plane symmetry (Figure 01 
we observe that for low temperature the vortex density 
augments when the tip passes over the surface. Initially 
the vortex density grows reaching quickly a constant av- 
erage. For higher temperature the vortex density is al- 
most insensitive to the tip indicating that the energy 
transfer becomes more difficult. At low temperature it is 
easier to excite the vortex mode since they have low cre- 
ation energy due to the out-of-plane spin component. At 
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FIG. 4: Vortex Density as a function of time for Dh = 
Ds = O.IJ and Xs = Xs ~ U- The graphics from left 
to right and top to bottom correspond to the temperatures 
T = 0.1, 0.2, 0.3, 0.4, 0.5, 0.6. The system has an out-of-plane 
symmetry. 



higher temperature the system is already saturated and 
creating a new excitation demands more energy. For the 
in-plane symmetry (Figure [SJ the situation is opposite. 
At low temperature there is no energy transfer to vortex 
modes. Creating a vortex demands much energy because 
the spin component of the vortex is almost whole in- 
plane. At higher temperature the system is soft so that it 
can absorb energy augmenting the vortex density. Even- 
tually it reaches saturation at high enough temperature. 
For the case when the vortex has an in-plane-symmetry 
(A < Ac), shown in figure|3 the average vortex density is 
constant even at higher temperatures. For (A > Ac) (Fig- 
ure [7|l the situation is similar to that where the system 
has a global out-of-plane symmetry (Compare to figure 
EI). 

In figure |S1 we show the velocity of the head and the 
vortex density as a function of time plotted in the same 
graphic for some interesting situations. The vortex den- 
sity is multiplied by a factor of 2 as a matter of clarity. 
In the first plot we observe that the kinetic energy of the 
tip is transferred to the surface. The head stops, moves 
back and forth and escapes from the surface influence. 
At higher temperatures the situation is a bit more com- 
plex. Depending on the initial condition the head passes 
through the surface region just augmenting the vortex 



FIG. 5: Vortex Density as a function of time for Dh = Ds — 
—O.IJ and Xs — Xs = U. Temperatures are the same as in 
figure 01 The system has an in-plane symmetry. 
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FIG. 6: Vortex Density as a function of time for X^ — Xs = 
0.6J and Dh = Ds = 0. Temperatures are the same as in 
figure 0] The vortex has an in-plane structure. 




FIG. 7: Vortex Density as a function of time for \s — Xs ~ 
0.9J and Dh = Ds = 0. Temperatures are the same as in 
figure 01 The vortex has an out-of-plane structure. 



density JHlb) or it can be trapped in the surface region, 
as seen in figure |S|c. The cost for the increase in the 
vortex density is a decrease in the kinetic energy. 

In figure El a and b we show the velocity and vortex 
density as a function of time for two different initial ve- 
locities {v = 0.6,0.2 respectively), at the same tempera- 
ture T — 0.2. For the higher velocity, the tip decreases 
its velocity almost up to stop, however, its kinetic energy 
is sufficient to go across the surface region. For a lower 
velocity, figure 05, the tip collides elastically with the 
surface. Because there is no lost in kinetic energy the 
vortex density is conserved. We note that the increase 
in the vortex density it is not an instantaneous response 
to the diminishing of the kinetic energy of the tip. It 
may be due to an intermediate mechanism: The kinetic 
energy is used to excite spin waves in the surface. Be- 
cause there is no mechanism of energy dissipation, part 
of the energy contained in the spin waves is transferred 
to vortex excitations. 



IV. CONCLUSIONS 

We have used Monte Carlo and spin dynamics simula- 
tion to study the interaction between two magnetic mo- 
bil parts: a magnetic reading head dislocating close to 
a magnetic surface. Our interest was to understand the 
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FIG. 8: Instantaneous velocity and vortex density as a func- 
tion of time for Dh = Ds = O.IJ and As = As = IJ. The 
full line is the vortex density and the dashed line is for the 
instantaneous velocity. From top to bottom T — 0.1, 0.2. The 
initial velocity is ii = 0.5 in all cases. The vortex density is 
multiplied by a factor 2 as a matter of clarity. 



mechanism of heat dissipation when the forces involved in 
the system are magnetic in essence. To simulate the sur- 
face we have considered a magnetic crystalline substrate 
interacting magnetically with a magnetic tip. From the 
results presented above we have strong evidences that 
vortices play an important role in the energy dissipation 
mechanism in magnetic surfaces. The augmenting of the 
vortex density excitations in the system increases its en- 
tropy. That phenomenon can blur any information even- 
tually stored in magnetic structures in the surface. An 
interesting result is the velocity behavior of the tip pass- 
ing close to the surface. In principle we should expect 
that the velocity will always diminish, as an effect of the 
interaction with the surface. However, for certain initial 
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conditions the effect is opposite. The tip can oscihate, 
be trapped over the surface or even be repeUed. In the 
case of an elastic cohision the vortex density remains un- 
changed. If the vortex density increases the tip's kinetic 
energy diminishes. However, the increase in the vortex 
density is not an instantaneous response to the diminish- 
ing of the kinetic energy of the tip. We suspect that an 
intermediate mechanism involving spin wave excitations 
is present intermediating the phenomenon. 

The effects on friction observed in our simulations 
demonstrate that when pure magnetic forces are involved 
they are quite different from ordinary friction. There are 
two points that should be interesting to study: The ef- 
fect of normal forces applied to the system and how the 
observed effects depend on the contacting area between 
surfaces. 



FIG. 9: Instantaneous velocity and vortex density as a func- 
tion of time for Dh = Ds = O.IJ and Xs = Xs ~ U- The 
full line is the vortex density and the dashed line is for the 
instantaneous velocity. The temperature is T = 0.2 in both 
cases. From top to bottom the initial velocity is i; = 0.6, 0.2. 
The vortex density is multiplied by a factor 2 as a matter of 
clarity. 
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